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Abstract 



In this paper, we study a vector conservation law that models the growth and selection 
of ovarian follicles. During each ovarian cycle, only a definite number of follicles ovulate, 
while the others undergo a degeneration process called atresia. This work is motivated by 
a multiscale mathematical model starting on the cellular scale, where ovulation or atresia 
result from a hormonally controlled selection process. A two-dimensional conservation 
law describes the age and maturity structuration of the follicular cell populations. The 
densities intersect through a coupled hyperbolic system between different follicles and 
cell phases, which results in a vector conservation law and coupling boundary conditions. 
The maturity velocity functions possess both a local and nonlocal character. We prove 
the existence and uniqueness of the weak solution to the Cauchy problem with bounded 
initial and boundary data. 
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1 Introduction 

In this paper, we study the following vector conservation law 
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$ft(t, x, y) + (A f cj)f(t, x, y)) x + (B f <j) f (t, x, y)) y = C(f f (t, x, y) 




te[0,T), (x, y) € [0, l] 2 



where $ f = ($[,■■■ ,~$n, <£{,••• A f N , 4>{,--- ,<fi f N ) T , / = 1> ■ ■ ■ ,n, and 

N N N 

A f := diag (g f , ■ ■ ■ ,g f , gf,--- ,?/, <7jV ■ ,<?/}> 

at 

5/ := diag {7i/, • • • ,7i/, 0, • • • ,0, /i f , • • • ,h f }, 
_JV_ jv N 

C := -diag { A, • • • ,A, 0~^~0, A, • • ■ ,A}. 

Here gf and are positive constants, ~g~f = gAuf) € C 1 ([0,oo)), hf = hf(y,Uf), 
hf = hf(y,Uf), A = X(y,U) and A = X(y,U) all belong to C 1 ([0, 1] x [0, oo)), with Uf = 
Uf(Mf(t), M(t),t) the local control and U = U(M(t),t) the global control, where Uf > and 
U > are continuous differentiable, i.e., u f G C^fO,©©) 2 x [0,T]) and ?7 € C^fO, oo) x [0,T]). 
For instance, the specific case that motivated our work is presented in Appendix 6.1. The 
function Mf is the maturity on the follicular scale given by 

N ,x ,i „ N .i r i 



M f( f ) := y2 / a i r 7sy4 ) k(t,x,y)dxdy + V / / (a 2 - ai)j^y4> f k (t,x,y) dx dy 
k=i Jo Jo k=i Jo Jo 

N ri ri 

+ Y1 / / a 2lo(loy + ls)(p{{t,x,y)dxdy 
k=1 Jo Jo 

and 

n 

M(t) :=X)M/(t) 
/=i 

is the global maturity on the ovarian scale. The parameters ai,a2,7s and 70 are positive 
constants with a 2 > a±. 

This work is motivated by problems of cell dynamics arising in the control of the develop- 
ment of ovarian follicles. Ovarian follicles are spheroidal structures sheltering the maturing 
oocyte. During each ovarian cycle, numerous follicles are in competition for their survival. 
However, only very few follicles reach an ovulatory size, most of them undergo a degenera- 
tion process, known as atresia (see for instance [E]). A mathematical model, proposed by 
F. Clement [5] using both multi-scale modeling and control theory concepts, describes the 
follicle selection process. For each follicle, the cell population dynamics is ruled by a con- 
servation law, which describes the changes in cell age and maturity. Two acting controls, 
Uf and U (see [T] and also Appendix 6.1), are distinguished. The global control U results 
from the ovarian feedback onto the pituitary gland and impacts the secretion of the follicle 
stimulating hormone (FSH). The feedback is responsible for reducing FSH release, leading to 
the degeneration of all but those follicles selected for ovulation. The local control Uf, spe- 
cific to each follicle, accounts for the modulation in FSH bioavailability related to follicular 
vascularization. In addition, the status of cells are characterized by three phases. Phase 1 



2 



and 2 correspond to the proliferation phases, and Phase 3 corresponds to the differentiation 
phase (see Fig 1). In Appendix 6.2, we have reformulated the original model to a new system 
(jl.ip . where the unknowns are all defined on the same domain [0,T] x [0, l] 2 , so that it can 
be treated as a general model for multiscale structured cell populations (see Appendix 6.2 
for the relation between the original notations and the new notations). 
The initial conditions are given by 



4>fo ■ =(ff(0,x,y) 



(1.2) 



--(4>io(x,y),- ■ -,4> m (x,y), (j){ (x,y),- ■ -,(j) S m (x,y), 4>{ (x,y),- ■ -,4> f m (x,y)f 



We use the simplified notations u f (t) := u f (M f (t),M(t),t), U(t) := U(M(t),t), u f (t) :-- 
u f (Mf(t),M(t),t) and U(t) := U(M(t),t) in the whole paper. 
The boundary conditions at x = are given by 



<p({t,0,y)=0, (t,y) G [0, T] x [0,1], 

^(M, !/),(<,!/) G [0,T] x [0,l],fc = 2 ) -- - ,N. 



4(t,0,y) = 

ai9_f(uf(t)) 

p k (t,o,v) = ai9f [ Uf{t)) 4 f k (t,i,y),(t,y) e [0,21 x [o,i],fc = ,N. 



(1.3) 



,N. 



T 9f 

4(t,0,v) = 0, (t,y) G [0,T] x [0,1], 
[ <t> f k {t,0,y) = <f> f k _ 1 (t,l,y), (t,y) G [0,T] x [0, 1], fc 

The boundary conditions at y = are given by 

4(t,x,0) = ${{t,x,0) =0, (i,x) G [0,T] x [0,1], fc = l, 



,N. 



d(t,-x,l), (t,x)€[0,T]x[0,^], 
ai 

0, (t, x) G [0, T] x [— , 1], 
a 2 



iV. 



;i.4) 



The boundary conditions at y = 1 are given by 

^(t,ac,l) = 0, (t, x) G [0, T] x [0, 1], 



k = !,■■■ ,N. 



(1.5) 



The well-posedness problems of the hyperbolic conservation laws have been widely studied 
for a long time. We refer to the works [H [T5l [T6] (and the references therein) in the 
content of weak solutions to systems in conservation laws, and [13\ Q3] in the content of 
classical solutions to general quasi-linear hyperbolic systems. In this paper, we perform the 
mathematical analysis of this system: we prove the existence, uniqueness and regularity of 
the weak solution to the Cauchy problem defined by (|l.ip - (|1.5p with initial and boundary 
data in L°°. The main difficulty tackled with in this paper comes from the nonlocal velocity, 
the coupling between boundary conditions and the coupling between different follicles in the 
model. Additionally, we have to deal with loss terms, and the problem is a 2 - space dimension 
one. 
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Figure 1: The illustration of N cell cycles for follicle /; x denotes the age velocity and y 
denotes the maturity velocity. The top of the domain corresponds to the differentiation 
phase and the bottom to the the proliferation phase 



In the related works that have considered nonlocal velocity problems [3J H] , the problems 
are only 1 - space dimension and do not have source terms. The velocities there are always 
positive, while in our case, the velocities hf(f = 1, • • ■ , n) we considered in this paper change 
sign in Phase 3. In another related work [10] . which was also motivated by [5j, the author has 
performed a mathematical analysis on this kind of model. He has reduced the model to a 1- 
space dimension mass-maturity dynamical system of coupled Ordinary Differential Equations, 
basing on the asymptotic properties of the original law. Once reduced, the model is amenable 
to analysis by bifurcation theory, that allows one to predict the issue of the selection for one 
specific follicle amongst the whole population. And he gave some assumptions on velocities 
according to biological observations. In another work the author also studied the well- 
poseness of the model, he proved the existence of a measure valued solution without proving 
the uniqueness of the solution, and he gave the behavior of the solution to the main equation 
along its characteristics. An associated reachability problem has been tackled in [6], which 
described the set of microscopic initial conditions leading to the macroscopic phenomenon of 
either ovulation or atresia in the framework of backwards reachable sets theory. The authors 
also performed some mathematical analysis on well-poseness of this model in a simplified open 
loop like case, there the authors assumed that the local control Uf and the global control U 
are given functions of time t. While in this paper, we consider the close/open loop case. 

This paper is organized as follows. In Section 2 we first give the main results. In Section 
3 after giving some basic notations, we prove that the vector maturity M := (Mi, • • • ,M n ) 
exists as a fixed point of a map from continuous function space, and then we construct a 
local weak solution to the Cauchy problem defined by (jl.ip - (|1.5p . In Section 4 we prove 
the uniqueness of the weak solution. Finally in Section 5 we prove the existence of a global 
solution. In complement, we introduce in Appendix 6.1 the original mathematical model 
proposed by F. Clement [5]. In Appendix 6.2 we reformulate this model to the new system 
(|1.1|) and we give the equivalence between the original notations and the new notations. In 
Appendix 6.3 we give a basic lemma that is used to prove the existence and uniqueness of 
the weak solution. 
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2 Main results 



We recall from [21 Section 2.1], the usual definition of a weak solution to the Cauchy 
problem defined by fjl.ip ~ fjl.5p . 

Definition 2.1. Let T > 0, <f> fQ € L°°((0, l) 2 ) be given. A weak solution of Cauchy problem 
(HU-flrSJ) is a vector function $ f G C°([0, T]; L 1 ((0, l) 2 )) such that for every r G [0,T] and 
every vector function (p :=(lp 1 ,- • ^i>' ' "> <?N ' "j <Pn) T £ -^([O, t] x [0, l] 2 ) with 

y?(r,x,y) = 0, V(x,y) € [0,1] 2 , (2.1) 

£(t,l,y) = 0, V(t,y)G [0,r] x [0,1], (2.2) 

¥>i(t,0,y) = £i(i,0,i/) = 0, V(t,y) £ [0,r] x [0,1], (2.3) 

? fc (t, x, 0) = fo(t, x, 0) = <^(i, x, 1) = £ fc (t, x, 1) = 0, V(t, x) G [0, r] x [0, 1], (2.4) 

one has 

mi 
H(t,x,y) ■ (0t(t,x,y)+Af<p x (t,x,y)+Bffi y (t,x,y)+C0(t,x,y))dxdydt (2.5) 

„i „i TV „ T „£1 



+ / / $fo(x,y) • <p(0,x,y)dxdy+y~] / / 2 h f (0,u f (t))(f){(t,—x,l)fk(t,x,0)dxdt 
Jo Jo Jo Jo a i 

+5_ / / ^ f k -i(tA,y) 1 PkitAy)dydt+y2 / / g f 4- 1 (t,l,y)^ k (t,0,y)dydt 

,._Jo J0 fl l 7 r, •/ J 



fc=2 " u " u x fc=2 

+ J2 r r 1 a 19fgf ( Uf (t)) - f hy) ~ k{tAy)dydt = 



With the definition, we have the main result 

Theorem 2.1. Let T > 0, <^>/o £ £°°((0, l) 2 ) &e given. Let us further assume that 

g~f(uf) > 0, Vii/G [0,oo), 
hf(y,u f )>0, V e [0,1] x [0,oo), 

hf(0,u f ) > 0, hf(l,Uf) < 0, Vwje [0, oo). 

T/ien i/ie Cauchy problem defined by U. 1\ )- [T75\) admits a unique weak solution <f>f = 
((/>{,••• , <j) N , </>{,••■ ,4>n, 4>{-,"' i4>n) T - Moreover, the weak solution <ff even belongs to 
C°([0,T];LP((0,1) 2 )) for allpe[l,<x>). 

The sketch of the proof of Theorem 12.11 consists in first proving that the maturity M(t) = 
(Mi(t),--- ,M n (t)) exists as a fixed point of the map M i— > G(M) (see Section 3.1), and 
then in constructing a (unique) local solution (see Section 3.2 and Section 4), before finally 
proving the existence of a global solution to the Cauchy problem defined by (jl.ip ~ fjl.5p (see 
Section 5). 
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3 Fixed point argument and construction of a local solution 
to the Cauchy problem 



In this section, we first derive the contraction mapping function G. Given the existence 
of fixed point to this contraction mapping function, we can then construct a local solution to 
the Cauchy problem defined by (jl.ip - (jl.5p . 

3.1 Fixed point argument 

First we introduce some notations. Let us define: 

ess sup (X:y)e[0A] 2\(f) k0 (x,y)\, 
ess sup {Xty)e[0A] 2\4> f k0 (x,y)\, 
ess sup iXty)e[0A] 2\4>{ Q (x,y)\. 

n N 

K 1=2^(70 + ls) 2 Y!Y^ ( a iPfco!lLi((o,i) 2 )+( a 2 - °i)llSolli, 1 ((o,i) a )'+ <l 2llSolU 1 ((o J i) 3 )) ) 

/=lfc=l 

(3.1) 

K x := max [g f , \\g f (u f (M f ,M, *))|| c i (gi) , \\h f (y, u f (M f , M, t))\\ c i (Q2) , (3.2) 

\\hf (y,u f (M f , M, t)) \\ci m , \\\(y, U(M, t)) \\ C i m , \\ X(y, U(M, t))\\ cHQa )}, 

K 2 :=mm\ inf g f (u f (Mf,M,t)), inf hf (y, ut (Mt, M, t))\ > 0, (3.3) 
I (Mf,M,t)eQi n 1 (y,Mf,M,i)eQ2 S M 1 J V ' 

Qi := [0,K] 2 x [0,T], Q 2 := [0,1] x [0,K] 2 x [0,T], Q 3 := [0,1] x [0, K) x [0,T]. 
For any 5 > 0, let 

ni,*:={#(*)=(Mi(t) ,.. ., Mn {t)) 6 C°([0, S\) : ||M|| c o (M) :=rnax \\M f \\ c o {m) <Kj , 

where the constant K is given by (|3.ip . 

In order to derive the expression of the contraction mapping function G, we solve the 
corresponding linear Cauchy problem f)l.l[) - f)1.5|) with given M £ ft>s,K- F° r an Y fixed t G [0, 5] 
and / € {1, • • • , n}, we trace back the density function <j)f at time £ along the characteristics 
to the initial and boundary data, hence we divide the plane time t into several parts. For 
Phase 1, the velocities hf are always positive, we introduce three subsets u;^'* an d wjj'* 
of [0, l] 2 (see Fig 2 (a)) 

cj{'* := |(x,y)| ^ g f {u f (a))da < x < 1, *7(*>^ 9 f(uf(cr)) da) < y < lj, 

wf* := {(a?, y)| 0<x< g f (u f (a))da, rj(t,x) < y < l}, 
E#* := [0,l] 2 \(znfu^). 



Il^fcoll : ~ H</>fcollL°°((0,l) 2 ) : ~ 

ll&ll := ll^fcolli-((o,i) 2 ) := 
fell := ll^olli-((o,i) 2 ) := 
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Here y = rj(t,x) (see Fig 2 (a))satisfies 

drj 



^L = h f (rj,u f )(s), rj(6) = 0, 



<s<t, 



(3.4) 



with 9 defined by x = / g t(uf(a))da. 

Je 

If (x, y) G U i , we trace back the density function <j) k at time t along the characteristics 
to the initial data. Otherwise, we trace back the density function at time t along the char- 
acteristics to the boundary data. For any fixed t G [0,(5] and (x,y) G [0, l] 2 , let us define 
characteristics £i( s ) := (xi(s),yi(s)),i = 1, • ■ ■ ,4 (see Fig 2), which will be used to construct 
the contraction mapping function. 

In the whole paper, we denote by (0,xo,yo) the points on the bottom face (t = 0), 
(to, ctQ, 0) the points on the left face (y = 0), (to, ao, 1) the points on the right face (y = 1), 
(to,0,/3o) the points on the back face (x = 0) and (to,1,A)) the points on the front face 
(x = l). 
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Figure 2: For Phase 1, time t plane is divided into three parts u>{'*, uJg an d wjj . (a) Case 
(t,x,y) G characteristic £2 connects (t,x,y) with (0, xo>2/o); ( D ) Case (t,x,y) G ZJ2 > 
characteristic £3 connects (i, x, y) with (to, 0, /3o); According to the coupled boundary between 
Phase 1 at x = 1 and Phase 2 at x = 0, we define characteristic £4 which connects (tq, 1, A)) 
with (0,x ,yo)- 

For any fixed x G [0, 1], we define characteristic £1 = (xi, y%) by (see Fig 2 (a)) 
^ =9f(uf(s)), ^ = hf(yx,u f )(s), £i(f) = (x,l). 



If (x,y) G cJ{' (see Fig 2 (a)), we define characteristic £2 = (^2, 2/2) by 



-£7 = 9f(uf(s)), -j^ = hf(y 2 ,u f )(s), &(t) = (x,y). 
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One has f 2 (s) e [0, l] 2 , Vs G [0,i]. Let us define 



(x ,y ) := (x 2 (0),y 2 (0)). (3.5) 
If (x,y) G uJ^' (see Fig 2 (b)), we define characteristic £3 = (£3,2/3) by 

"J^ = 9f(u f (s)), -j^ = h f (y 3 ,u f )(s), &(t) = (x,y). 
There exists a unique To € [0,i] such that £3(70) = 0, so that we can define 

A) := Jtefo). (3-6) 

According to the coupled boundary between Phase 1 at x = 1 and Phase 2 at x = 0, for any 
fixed (xo,yo) ^ [0>1] 2 > we define characteristic £4=(x4,y4) by (see Fig 2 (b)) 

=9f(uf(s)), ^ = hf(y 4 ,u f )(s), £ 4 (0) = (x ,y ). 

For Phase 3, due to the fact that velocities fry change sign in Phase 3, we divide the plane 
at time t into four subsets tjf'*, wjj , SJ3 and a3{'* of [0, l] 2 (see Fig 3 (a)) 

: =|(x,y)|5/t < x < 1, rjx(t,g f t) < y < rj 2 {t,gft)}, 

:={(x,y)|0 <x<g f t, < y < ^i(t,x)}u{(x,y)| ^ < x < 1,0 < y < ^(t.flfft)}, 
^3'* : ={( s >y)l < z < r/i(i,x) < y < %(i,^)}, 
:=[0,l] 2 \(5fu5 2 / ' t U^' i ). 

Here y = rji(t,x) and y = r)2{t,x) (see Fig 3 (a)) satisfy 

dfj^ ~ x x 

-j- = hf(rji,Uf)(s), Vi(t-—) = 0, t- — <s<t, 
ds 9f 9f 

dm 7 ,~ \ / \ ~ , x . x 

— = hf{r)2,Uf)(s), n 2 {t-—) = l, t- — < s <t. 

ds 9f 9f 

If (x, y) G w^' , we trace back the density function <p{ at time i along the characteristics to the 
initial data. Otherwise, we trace back the density function at time t along the characteristics 
to the boundary data. For any fixed t G [0, 5] and (x, y) G [0, l] 2 , let us define characteristics 
£i(s) := (xi(s), yi(s)),i = 5, • • • ,10 (see Fig 3), which will be used to construct the contraction 
mapping function. 

If (x,y) G uj[' (see Fig 3 (a)), we define characteristic £5 = (xs,ys) by 

~j^ = 9fi = hf(yB,Uf)(s), = (x,y)- 

One has £ 5 (s) G [0, l] 2 , Vs G [0,*]. Let us define 

(x ,ift)):=(s 5 (0), 1/5(0)). (3.7) 
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If (x,y) G u>2 (see Fig 3 (a)), we define characteristic ^6 = (xe,ye) by 
~J~ = 0/' ^ =hf(y6,u f )(s), €e(t) = (x,y). 
There exists a unique to G [0,t] such that j/e(*o) = 0, so that we can define 

"o := x 6 (t ). (3.8) 

For any fixed cell cycle k = 1, • • • , N, according to the coupled boundary between Phase 
1 and Phase 3 (corresponding to (p k (t,x,y) at y = 1 and cf) J k (t,x,y) at y = 0), we separate 
the right face of Phase 1 into two parts. For k = 1, • • • ,N, we define characteristic £7 passing 
through the right face intersects with the bottom face at (0, xo,yo) (see Fig 3 (b)). For 
k = 2, • • • ,N, we define characteristic £s passing through the right face intersects with the 
back face, and then back to the bottom face at (0, xo,yo) of the k — 1 cell cycle in Phase 2 
(see Fig 3 (b)). Hence, for any fixed (xo,yo) G [0, l] 2 , we define £7 = (£7,1/7) by 

-J^ =9f(uf(s)), ^ = hf(y 7 ,Uf)(s), £7(0) = (x ,y ), 

and £ 8 = (xs,y&) by 

^7 =9f(uf(s)), ^ = hf(yn,Uf)(s), ^ ^ X ° ) = (0,y ). 




Figure 3: For Phase 3, time t plane is divided into four parts o;{ ,£ , , £^3 and . (a) 

r ^ ft 

Case (t,x,y) G u^' , characteristic £5 connects (t,x,y) with (0,xo,yo); Case (t,x,y) G u;^' , 
characteristic £q connects (t,x,y) with (to,ao,0); (b) According to the coupled boundary 
between Phase 1 at y = 1 and Phase 3 at y = 0, we separate the right face into two parts, 

■~— f t 

and we define characteristic £7 and characteristic £§; (c) Case (t,x,y) G oj s ' , characteristic £9 
connects (t, x, y) with (to, 0, /3o); According to the coupled boundary between two consecutive 
cell cycles, we define characteristic £10 which connects (to, 1, /3q) with (0, xo,yo)- 
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If (x,y) € (see Fig 3 (c)), we define characteristic £9 = (xg,yg) by 
dxg „ dyg ~ 

~fo = 9f> ~d7 = h f(y9,uf){s), £g[t) = (x,y). 

There exists a unique To € [0, t] such that xg(ro) = 0, so that we can define 

A) := Ifefo). (3.9) 

According to the coupled boundary between two consecutive cell cycles (corresponding to 
0jfc-i(*j x, y) at x = 1 and <^(i, x, y) at x = 0) (fc = 2, • • • , iV), for any fixed (xq, yo) € [0, l] 2 , 
we define £i = (xio,yio) by (see Fig 3 (c)) 

~jT = 9fi ~xf = hf(yio,Uf)(s), £10 (0) = (x ,y ). 

With all of the above defined characteristics and noting Lemma 16. II in Appendix 6.3, we 
are now able to define a map G(M)(t) := \ Gi(Mi,M)(t), ■■■ , G n (M n , M)(t)j for all t € [0, 5} 
with 

pi(0) (■l-J*g f (u f (a-))da N t _ 

G f (M f , M)(t):= / / a l7 2 £ A(*o, VoMfy- & X ^' U){ - s) ds dx Q dy (3.10) 

Jo Jo k=1 

• l mC 1 ^) N _ -!\- X0 \{y 3 ,U){s)ds 



+ / / 2(a 2 - ai)7 2 (x , y )y 3 (*)e ^ dy dx 

Ji-gftJo k=2 

rl rl-9ft N _ 

+ / / (o 2 - aih 2 Y\ </>£ (x , Vo)dx dy 
Jo Jo k=1 

+ / / a l7 2 V(yo+/ ^(y4,«/)(a)da)^ (xo,yo)e- i ° T0 ^ 4 ' [/)(s)d %o^o 

Jl-f*g f (uf(a))da-J0 , , 7o 



1 



+ 



JO 



AT 

a 2 7o^(7oy5(t)^)^ ^o,yo)e- / o^fe' c/ )( s ) d ^x d 2 /o 



fe=i 



+ f /'"ana E(7oye(t)^)^ (xo, yo)^ ^ , yo 

Jyiioyo k=1 

rl rl N ^ -iS[°- xn KV8,U)(s) ds+Jl X(y 6 ,U)(s) ds) 

+ / / , 2(a 2 -ai)7o V"(7oy6(i)+7 s )^_ 10 (x ,yo)e 9 -f dy dx 

1 ~ 3: 

rl rl N-l _ -(f s f X( yi0 ,U)( s )ds+ji- Xn \(y g ,U)(s)ds) 

+ / / a 2 7o y^(7oy9(i)-^s)^ ( x o,yo)e 9 ? dx dy . 

J0 Jl-gft k=l 

Here y%{t) (i = 2, • • • ,10) can also be determined by (xo,yo), to is defined by 1 — xo = 

/•TO 

/ ~g~f(uf(a))da, to is defined by yj(tg) = 1, and to is defined by ys(to) = 1- In <|3.10|) . 
Jo 

yo = C( x o) '■= ^(Oj^o) (see Fig 3 (b)), where r](s,xo) satisfies 

^ = hf(v,Uf)(s), v(0) = h 0<s<9, 
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f _ 

with 9 defined by 1 — x = / g j(uf(a))da. 

Jo 

Next we prove the following fixed point theorem. 

Lemma 3.1. If 5 is small enough, G is a contraction mapping on £lg,K with respect to the 
norm. 

Proof: It is easy to check that G maps into ^l^,K itself if 

0<<5<min{JL T}, (3.11) 

where K x is defined by (|5^|> . Let M(t)=(Mi(t),- ■ -,M n (t)), M {t)={M x {t) • -,M n (t)) GO^, 
and M=Y^} =1 M f , M=£" =1 M/. In order to estimate \\G(M) - G(M)|| c o( [0 we first 
estimate the norms \\Gf(Mf,M) — Gf(Mf,M)\\ C! on m separately Observing the definition 
(|3,10p of Gf, it is sufficient to estimate — ?/j||c ([o,<5])> where 

rPi 

Vi = C+ / hf(yi,u f )(a)da, 



Vi = C+ h f {y^u f ){a)da. (3.12) 

J OLi 

Here C denotes various constants, ht represents ht or hf, and 0<ai<Pi<t<5<mm{ , T}. 

2Ki 

We have 

\Vi(s) - yi(s)\ < / \h f (yi,u f )((T) - h f (yi,u f )(a)\da 

J di 

< tKx{\\Mf - M/ 1 1 c?o ([0,5]) + W - M\\ c o {m) ) + tKxWyt - yi\\c°{[ ai M> 

licncc 

tiri(||M / -M / || c o ([0i5]) + ||M-M|| c o ([0 . 5]) ) 

Ik - i/i||co([o,«]) < i~tK ' ^ ' 

By (|3.13p . we have 



\hf(yi,Uf)(a) - hf(yi,u f )(a)\da (3.14) 

< tifi(||M/ - M f \\ c u im) + \\M - M|| c o ([0)5 ])) + - 2/i|lc°([o,<5]) 

<tK(\\M Mil +IIM Mil ^fC^M / || c o ([0 , 5]) +||M-M|| a o (M) ) 



1 - tKi 

_ mjWMf - M f \\ c0[m) + \\M - M\\ c o {m) ) 
1-tKx 

Similarly we have 

\X(Vi,U)(a) - \{y h U){a)\da <tK x \\M - M\\ c o {m) + tKx^ - yi\\c°([o,S]) 

t 2 Kf\\M f - M f \\ c0{m) + mWM - M|| c o (M) 



< 



1 - tKi 

(3.15) 
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where A represents A or A, and 

\9f(u f (a))-g f ( Uf (a))\da < ^i(||M/ - M f \\ c0(m) + \\M - M\\ c o {m) ). (3.16) 
By (j3~T3]) - fl3~T6D and the definition (fXTU]) of Gf, we have 
\\Gf(M f ,M)-G f (M f ,M)\\ CO{m) 

< tC[\\M f - M f \\ c o (m) + tC[\\M - M\\ c o {m) . (3.17) 
The expressions of C{ and C{ are given by (|6.30p in Appendix 6.3. They depend on 

Pfcoll> ll^feoll' ll^fcoll'-^l and K %- TllUS We have 

||G(jtf)-£(M)|| c o (M ) (3.18) 
= (g 1 (M 1 ,M)-G 1 (M 1 ,M),--- ,G n (M n ,M)-G n (M n ,M)) 

\ / C°([0,o]) 

= max \\G f (M f ,M)-G f (M f ,M)\\ c o [m) 



< max 



Lax (tC{\\M f - M f \\ c0[m) +tC f 2 \\M - M|| c o (M) ) 



< max [tC[ \\M f -M f \\ CO{m) +tC> \\Mi-M 1 \\ c o {M) +- ■ ■ +tC } 2 \\M n -M n \\ c o {[m 
<maxt(C{ + C Y 2 / )(||Mi - Mi|| c o ([0)(5]) + • • • + \\M n - M n \\ c o {m) ) 
<ntmax(C{ + C|) max \\M f - M f \\ c o {m) . 



We finally get 

\\G(M) - G(M)\\ c o {m) < ntmax(C{ + C()\\M - M\\ c0{m y (3.19) 

f — ¥ ' — ¥ 

Hence we can choose S small enough (depending on ||0fcoll> ll^fcoll' H^i-oll' ^i, so ^at 

\\G(M)-G(M)\\ CO{m) < - M\\ CO(m) . (3.20) 

□ 

By Lemma 13.11 and the contraction mapping principle, there exists a unique fixed point 
M = G(M) in Q S>K . 

3.2 Construction of a local solution to the Cauchy problem 

Now we show how the fixed point M allows to find a solution to Cauchy problem (jl.lj) - 

(jl.5p for t € [0,6]. Let us recall the definition of three subsets tj{' 4 , tJ^'* and tJg'* of [0, l] 2 , 

the definition of the characteristics £2>£3 and also the definition ()3.5[) of (xo,yo) and (|3.6p of 

r 

(ro,/?o)- For k = 1, we define ^(t^x^y) by 

(P{(t,x,y) := { ^w{x ,y )e , it [x, y) E u) x , ^ 21 ^ 

0, else. 
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f 

For k = 2, • • • , N, we define 4> k (t, x, y) by 



j 



^ (x ,y )e- / « [I(s,2 ' C/)+ ^ (?/2 ^ )]{CT),iCT , if (x,y) G sf*, 



2t„ 



a l9f(Uf(T )) 

0, else. 



8a 



tier 



, if(x,y)G^, (3-22) 



Since the dynamic in Phase 2 amounts to pure transport equations, we define (f)l(t,x,y), 
k = !,■■■ ,N by 



<i>lo( x -9ft,y), if (z,y) e [<?/t,l] x [o, l], 

a\g~f{t — £-) . x 

^- ^(t__,l,y), if (x,y) G [0,5/t] x [0,1]. 

T 9f 9f 



(3.23) 



Let us recall the definition of the four subsets w^'*, and of [0, l] 2 , the definition 
of the characteristics £5, £6)^9 and also the definition (|3.7p of (xoj?/o)> (13.81) of (to,ao) and 
(|3.9p of (ro,/3o)- For k = 1, we define (j){(t,x,y) by 



^o(xo,y )e- / ''^ 5 ' C/)+ ^ (3/5 ^ )1(ff)dCT , if (x,y) G Sf*, 



4>{(t,x,y) := < 



^(to, — ao,l)e Jt o l 9 » ' ;JW , if (a, y) G , 

0, else. 



(3.24) 



For k = 2, • • ■ ,N, we define <?!>{(£, x, y) by 



4>{(t,x,y) := < 



' ^ (x ,y )e- / o [ " fe ' C/)+ ^ ( ^ )1(<T)dCT , if (x,y) G 
^(to,^«o,l) e -^^^) + *^^W^, if (x, y) G Sg'*, 

0{_ 1 (r o ,l,/3o)e"A^^ + ^^^)]M^, if (*,„) e 



(3.25) 



0, else. 



Next we prove that the vector function <j)f defined by (|3.2ip - (|3.25p is a weak solution to 
Cauchy problem (ll.ip - (ll.5p for t G [0, 5]. To that end, we first prove that 4>f defined by Q3.21|) - 
([335]) satisfies equality {23]) of Definition EH then we prove that $ f G C°([0, 5]; L x ((0, l) 2 )). 
Let r G [0, 5]. For any vector function (p G C 1 ([0, t] x [0, l] 2 ), (p := (Tp 1 , • • • , yjjy, £>i, • • • , (fN, 
<?!,-•• , <?Af) T such that dTTJ-dM]) hold, by definition (133T]) - (|3~25|) of 4> f , we have 



ml ' 
cj>f(t, x, y) ■ (<p t + Af(p x + Bjify + C<p)dxdydt := V] A*. (3.26) 
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In 

N 



i v /* T" f f dh 



fc=i 



W 2 ^ ^ o5l ^ o ^ ^ 

Ax=2_jJJ a - ^ ^ e J V ^55 W- ^ ((p kt +g flpkx +hf<p ky -\ip k )dxdydt 



-/'■' 



AT rr rl <ll 



.>,X/ /•/ 



9ft 



k=l 

N 



A 4 :=V7 / / ■ ^— cf> k (t- —,l,y)(ip k t+gflpkx)dxdydt, 



^=1 ~/,t 



fc=r 

AT 



z* 7 " r r t ~ dhf 

Let us consider the first term A\ as an instance. By the change of variable (x,y) — > (xo,yo) 
and noting (|6.21|) of Lemma 16. II in Appendix 6.3, we have 

N rr r /3 pa 



k=1 J0 jo jo K 

+ h f Tp ky (t, x 2 (t),y 2 (t)) - XTp k (t, x 2 (t),y 2 (t))^j dx dy dt, 

where (see Fig 4 (a) and (b)) 

a:=l- [ g f (u f (a))da := h(t), /3 := 1 - / h f ( yi , u f )(a)da := g x (t). (3.27) 
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Clearly, a is a function of j3, suppose that a = h{(3). After changing the order of integration, 
A\ can be rewritten as 

^=E{/ //+/// +/ / / +/ / / 

fc=1 V0 JO JO JO Jfi(r)J0 J gi (r)J0 J J gi(r)J h{/3)J 

,1 rl rg^ifi d (V fo^ U ^Tp k (t, X 2 (t), y 2 (t)j) 

+ 11 / \<t>koi x o,yo) — "77 L dtdx y 

Jf 1 (T)Jh- 1 (a)Jo ' at 

r9i(r)rfi(x)_f ra^ri _j 

=zJr / 4>ko( x o,yo)Vk( G i x o>yo)dx dyo- / (f) k0 (xo,y o )7p k (0,x ,y o )dx o dy 

t^i 1 JO JO JO Jf x (T) 

1 rh(r) _ rl rl _ 

4>ko( x o,yo) l Pk(0,xo,yo)dxodyo- / (t) ko (x ,y )(p k (0,x ,yo)dx dy 



gi{rV0 J gi {ryh{0) 

1 rl -f 

<t>ko( x o, yo)^fc(0, x ,y )dy dx 



pi(r) /•! 



+ 



/ — 1 ( ) — 

^Ux ,y )e-Io 1 " ^(^^)(-)^(/ r 1 (a), a ; 2 (/ 1 - 1 (a)),y2(/r 1 («)))^o^o 



J/i(r) 



+ / / ^o(^o,yo)e-^ 1 Afe2 ' c/)(CT) ^ fc (5r 1 (/3),x 2 ( ff r 1 (/3)),y 2 (<7r 1 (/3)))^od 2/ o 
+ / / ^ (x ,yo)e-i'o 1 " J ^ u ^Tp k (f{\a),x 2 (f{\a)),y 2 (f{\a)))dx dy 



+ 



4 (x ,y )e-^ 1 < " )A ^'^(^ fc ( 9 r 1 (/3),X2(5r 1 (/3)),y 2 (9r 1 (/3)))^o^o}. 

/i(r>//i-i(Q) J 



Changing the order of integration again (see Fig 4 (c)), we obtain 

1 r l _ rl rh(fl)_. _ 

4'ko( x o,yo) l Pk( (j ^ x o,yo)dyodx =- / (f) kO (x ,yo)(p k (0,x o ,yo)dxodyo. 



Mryh-^a) Jgi(r)JMr) 

(3.28) 



By (13381) and noting that x 2 (f{ 1 (a)) = y 2 {g^ l {0)) = 1, we get 

A i = ~^2 / ( Pko( x o,yo)^k( , x o,yo)dx dyo. (3.29) 



k=i 



o Jo 
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Similar to A\, we can prove that 



A, 



A* 



A, 



Ay 



N , T j 



E 

k=2 




It, 



if 2f 



o Jo °i 



0Li( T o> !> Po)<Pk( T o> °> Po)df3 dT , 



N ,1 ,1 



E 

fe=i 



JO 



^fco( x o, yo)^fc(0, x ,y )dx dy , 



N rT x 

E 



fe=i 




/( o)) ^(7-0,1^0)^(70,0,^0)^/80^0, 



j 



T 9/ 



^ ,1 ,1 



E 

fe=i 



iv 




J 



^fco( x O' 2/0)^(0, xo, ya)dx dy , 



' hf(0,u f (t ))(j) f k (to, —a , l)ip k (to,a ,0)da dt , 

JO a 2 



E 

fc=i 

AT „ T „x 

/ / 5/^{_i(To,l,^o)^fc(To,0,/3o)d/3odro. 



fc=2 



JO 



(3.30) 
(3.31) 
(3.32) 
(3.33) 
(3.34) 
(3.35) 



}(t,x,y) 



.(0.)!) \i 



1 1°. Wo) V 




(a) 



(b) 



(c) 



Figure 4: (a) For any fixed £ £ [0, r], the characteristics ^1,^2 and the definition of (b) 
The illustration of functions a = fi(t) and /3 = gi(t), t G [0, r]; (c) For changing the order of 
integration. 



By (|3.29|) - (|3.35p . we have proven that the vector function <j>f defined by (|3.21|) - (|3.25p 
satisfies (|2.5p . Next we prove that <fif G C°([0, 5]; L 1 ((0, l) 2 ))- Moreover, we can prove that 
0/ even belongs to C°([0, 5]; D>((0, l) 2 )), for all p G [1, 00). 

Lemma 3.2. The weak solution cpf to Cauchy problem (|l.ip - (|1.5p belongs to C°([0, 5]; £ p ((0, l) 2 )) 
/or aZZ p G [1, 00). 

Proof: From the definition (l3T2TT) - (l3^5D of 0/, we get easily that $ f G L°°((0,<5) x (0, l) 2 ). 
Next we prove that the vector function cfif belongs to C ([0, 5]; L p ((0, l) 2 )) for all p G [1, 00), 
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i.e., for every t,t G [0, 5} with t > t (the case that i < t can be treated similarly), we need to 
prove 

1 1 </>/(*> •) - -)IIlp((o,i) 2 ) °; as I* — *1 — ^ 0, Vjpe[l,oo). 

In order to do that, we estimate \\<f> k (t, •) - 4> k (t, ■) llz,p((o,i) 2 ): ~ llz,p((o,i) 2 ) and 

-)IUp((o,i) 2 )(^ = V - ,N) separately. 
Suppose that the characteristic passing through (t, 0, 0) intersects time t plane at (t, q±, q 2 ) 
(see Fig 5 (a)), we have 



<?i = / g f (u f (a))da, q 2 = h f (y,u f )(a)da, 
Jt Jt 

where y(s),t < s < t satisfies 

^ = hf(y,u f )(s), i/(t) = 0. 

We get easily that 

qi<C\t-t\, q 2 <C\t-t\. (3.36) 

Here and hereafter in this section, we denote by C various constants which do not depend 
on t, t, x or y. Let S := [q\, 1] x [q 2 , 1] (see Fig 5 (a)), following our method to construct the 
solution 4> k , when (x, y) G S, we have 

\4>l(t,x,y) -4>{(t,x,y)\ 

= ft<t,z,jfte-^n+%f9>>»M*)*' -4(t,x,y) 

< \4>{(t, x, y) - 4(t, x, y) \e~ J^^+^G&wMM *> + c$l{t,z,y)\\i-t\. 

Here £2 = (x 2 ,y 2 ) denotes the characteristic passing through (t,x,y) that intersects time t 
plane at (t,x,y). Hence (x,y) is defined by 



x = x 



9f(u f {a))da, y = y - J h f (y 2 ,u f )(a)da. 



We have 

\x-x\ < C\t-t\, \y-y\ < C\t-t\. (3.37) 

£ 

Since cj> k G L°°((0,6) x (0, l) 2 ), for every I > 0, there exists C; such that for every t G [0,6], 
there exists 4> k (t, •) G C 1 ([0, l] 2 ) satisfying 

\\4>h fa ') - 0fc(<> OlUf ((0,1)2) < y, ||(/> fc (t, -)llci([o,i]2) < Q. (3.38) 



Here and hereafter in this section, we denote by C\ various constants that may depend on I 
(the index of the corresponding e 
are independent of0<i<i<<5. 



(the index of the corresponding approximating sequences <p k (t, •), 4> k (t, ■) and <p k (t, •)) but 



17 



Thus, we have 



\4>k(t,%,y) ~ 0it(*>»>y)| 
<C\4> k \t,x,y) - ~4> k {t,x,y)\ + C\^ k (t,x,y) -~4> k (t,x,y)\ 

+ C\4 f k (t, x, y)\\t -t\ + Q\t - t\. (3.39) 



By (13381) and (13T39D . we have 



\<pl(t,x,y)-4(t,x,y)\Pdxdy<^ + C l \t-t\P. (3.40) 



Noting (|3.36p . we have 



\d(t,x,y) - <p{(t,x,y)\ p dxdy < C{ qi +q 2 ) < C\t - t\. (3.41) 



[0,i] 2 \S 

Combining (I3.40|) with (I3.41h . we obtain 



1 \4(t,^y)-d(t,^y)\ P dxdy<^- + C l \t-t\P + C\t-t\. (3.42) 




o jo 



IP 






(a) 



(b) 



(c) 



Figure 5: (a) For Phase 1, the definition of 5 and characteristic £2 which connects (t,x,y) 
with (t,x,y), characteristic passing through (t, 0, 0) intersects time t plane at (t, 91,92); (b) 
For Phase 2, the definition of S; (c) For Phase 3, the definition of S and characteristic £5 
which connects (t,x,y) with (t,x,y), characteristic passing through (t, 0, 0) intersects time t 
plane at (t, v, v\) and characteristic passing through (i, 0, 1) intersects time t plane at (t, v, ^2). 



Next we estimate 



4> k {t, •) ||iP((o,i) 2 )- F° r Phase 2, the characteristic passing 
through (t,0, 0) intersects time i plane at (t,gf(t — t),0). Let 5 1 := [(//(£ — t),l] x [0,1] (see 
Fig 5 (b)), when (x,y) £ S, we have 



h l(t,x,y) - (j) J k {t,x,y)\=\(l) J k {t,x,y) - <l> J k (t,x,y)\=\<f) J k (t,x -g f (t - t),y) - 4> J k (t,x,y)\ 



li- 



ft 
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Since 4> k G L°°((0,5) x (0, l) 2 ), for every I > 0, there exists C\ such that for every t G [0,6], 
there exists 4>{\t, •) G C 1 ([0, l] 2 ) satisfying 

OHl^o,!) 2 ) < y» H^(V)llcM[o,i] 2 ) < Cj- ( 3 - 43 ) 
Similar to (|3.40p . we can prove that 

II \^ k (t,x,y)-^ k (t,x,y)\Pdxdy<^ + C l \t-tr. (3.44) 
It is easy to check that 

\ft k (i,x,y) - ft k (t,x,y)\ p dxdy < C\t - t\. (3.45) 



[o,i] 2 \s 

Combining (|3.44|) with (|3.45j) . we get 



J (t, x, y) - ££(t, x, y)| p dxdy < ^ + C,|t- i| p + C|t - t|. (3.46) 



o jo 



Finally, we estimate \\4>{(t, •) — ')lli>((o,i) 2 )- Suppose that the characteristic passing 
through (t, 0,0) intersects time i plane at (t, v, V\), and the characteristic passing through 
(t, 0, 1) intersects time t plane at (t, v, v 2 ) (see Fig 5 (c)). Similar to (|3.36p . we can prove that 

v<C\i-t\, vi<C\t-t\, l-v 2 <C\t-t\. (3.47) 

Let S := [v, 1] x [-ui,^] (see Fig 5 (c)), when (x,y) G S, we have 

\4>{{l x, y) - g(t, x, y) I = \4>{{t, x, y)e~ ^^ + ^' n ^ a ^ -^ k {t,x,y)\. 

Here 2/5 = (£5,2/5) denotes the characteristic passing through (t,x,y) that intersects time t 
plane at (t,x,y). Similar to (|3.37[) . we can prove that 

\x - x\ < C\t - t\, \y - y\ < C\t - t\. 

Since 4>{ G L°°((0,5) x (0, l) 2 ), for every I > 0, there exists C\ such that for every t G [0,(5], 
there exists 0{ f (t, •) G C 1 ([0, l] 2 ) satisfying 

U{ l (t,-)-ft k (t,-)\\ LPmm < j, Ul\t,-)\\ CH[0 , 1]2) < Q. (3.48) 
Similar to ()3.40p . we can prove that 

II \^ k (t,x,y)-^ k (t,x,y)\ p dxdy<^ + C l \t-t\P. (3.49) 



S 



Noting (|3.47p . we have 

4>{(t, x, y) - %(t, x, y)\ p dxdy < C(v + v x + 1 - v 2 ) < C\t - t\. (3.50) 



[0,1] 2 \S 
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Combining (|3.49j) with (|3.50|) . we get 

J ${(t, x, y) - ft k (t, x, y)\ p dxdy <^ + Q\t- t\' + C\t - t\. (3.51) 

Therefore letting I be large enough and then \t — 1\ be small enough, the right hand side of 
(|3.42p . (|3.46p and (|3.5ip can be arbitrarily small. Above all, we prove that the vector function 
$ f belongs to C°([0, 5};L P ({0, l) 2 )) for all p G [1, oo). □ 

Let us recall Definition [2J] of a weak solution, we prove that the vector function <f>f defined 
by (|3.2ip - (|3.25p is indeed a weak solution to Cauchy problem (|l.ip - (jl.5p . 



4 Uniqueness of the solution 

In this section, we prove the uniqueness of the solution. 

Lemma 4.1. The weak solution to Cauchy problem (jl.ip - (jl.5p is unique. 

— =f =f —f —f ~f —f 
Proof: Assume that 4> f = {<j> x ,--- An^Ii' ' ' AnAi,--- , 4>n) T ^ C°([0, 5]; L 1 ((0, l) 2 )) is a 

weak solution to Cauchy problem (|l.ip - (|1.5p . Similar to Lemma 2.2 in [3l Section 2.2], we 

can prove that for any fixed t G [0, 5] and any 3>(t, x, y) : = • • • , $at, 

$i, • • • , 8jv, ii, • • • , § N ) T g c^qo, t] x [o, i] 2 ) with 

%,l,y)=0, V(r,y) € [0,t] x [0,1], 

$i (r, 0, y) = £ i (r, 0, y) = 0, V(r, y) G [0, i] x [0, 1] , 

¥ fc (r,ac,0) = * fc (r I ac,0) = * fc (r,s,l) = 8fc(r,x,l) = 0, V(r,x) G [0,i] x [0,1], 
one has 

m0/(r, x, y) • ($ r (r, x, y)+A/C(r, x, y)+B/ $ y (r, x, y)+C$(r, x, y))dxdydr (4.1) 
,i r i N rt p Bi 



+ / / <f f0 (x,y) ■ $(0,x,y)dxdy+y] / / " 2 /i/(0, «/(r))^ fe (r, — x, l)$ fc (t, x, 0)cfedr 

JO JO JO JO °1 

+ Y1 / / — h-i(r,l,y)^k{r,0,y)dydT+y2 / / g f cp k _ 1 (T,l,y)^ k (T,0,y)dydT 
k=2 JoJo °i k=2 JoJo 

+ aig/g ^ /(T)) ^(r,l,y)6 fc (r,0,y)dydr = f [' $ f (t,x,y) ■ $(t,x,y). 



o Jo 



In (|4.1|) . the velocity matrices Af,Bf and C are defined the same way as ^4/, i?/ and (7 but 
with Mf and M instead of Mf and M, where 



N f 1 f 1 =f N f 1 f 1 

M f( t )=y2 / a ilsV ( Pk{ t ^ x ^y)dxdy + / (a 2 -oi 

^ f 1 f 1 — / 

+ / / a 27o(7oy +^/ s )^ k (t,x,y)dxdy, 

1 ■/ ■/ 



hsV ( Pk( t ^x,y)dx dy 



20 



and M(t)=Y?} =1 Mf(t). Let <p (x, y)=((p 01 , ■ ■ ■ , p 0N , <Pox, ■ ■ ■ , ^ojv, &01, ■ ■ ■ , <Auv) T e C^((0, l) 2 ), 
and choose the test function as the solution to the following backward linear Cauchy problem 

" l r + Af$ x + B f $ y = -C$, 0<r<t, (x, y) G [0, l] 2 , 
£(i,x,y) = (fo(x,y), (x, y) G [0, l] 2 , 

%,l,y) = 0, V(t, y) G [0, t] x [0, 1], (4.2) 
¥i (r, 0, y) = $i (r, 0, y) = 0, V(r, y) G [0, t] x [0, 1], 

[ $ k (T,X,0)=$k(T,X,Q)=$ k (T,X,l)=$ k (T,X,l)=0,V(T,x) G [0, t] X [0,1]. 

For any fixed t G [0, <J], we introduce three new subsets cJ 2 '* an d of [0, l] 2 



:= |(x,y)| y g f (u f (a))da < x < 1, 9 f{uf{<?)) da) < y < l|, 

:= {(x,y)| < x < / g f (u f (a))da, fj(t,x) < y < l}, 



cj 2 

:= [0,1] 2 \(^ < U^ < ). 
Here y = r/(i, x) satisfies 

■£ = hf%u f )(s), W) = 0, 9<s<t, 
ft 

with defined by x = / gf(Uf(a))da. 

For any fixed t G [0,(5] and (x,y) G [0, l] 2 . If (x,y) G ujf'*, we define £ 2 = 0^2, 1/2) by 

=9f(u f (s)), ^ = h f (y 2 ,u f )(s), | 2 (t) = (x, y). 

Let us then define 

(x ,y ) := (x 2 (0),y 2 (0)). (4.3) 
If (x,y) G tJ 2 , we define £ 3 = (x 3 ,y 3 ) by 

=9f(u f (s)), ^ = h f (y 3 ,u f )(s), £ 3 (f) = (x, y). 

There exists a unique To such that X3(r"o) = 0, so that we can define 

/So-lfefo). (4.4) 

By ([HQ and (B~2j) . we obtain that 

n$ f (t, x, y) ■ <p f Q (x, y)dxdy (4.5) 
,1 ,1 JV 



= / / $ f0 {x,y) ■$(0,x,y)dxdy+J2[ /° 2 /i/(0, «/(t))^(t, — x, l)$ fc (t, x, 0)dxdr 
+Z7 / — <f>k-i(T,l,y)$k(T,0,y)dydT+J2 9f^k-i(rA,y)^k(r,0,y)dydT 
E/' /' - g/?/( ° /(T)) a!(r,l.,)8.(r,Q,l>)^r. 



N 

+ 
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Let us recall the definition of the characteristic £ 2 an d also the definition (|4.3p of (xQ,y ). 
By solving the backward linear Cauchy problem (|4.2[) and noting (|6.2ip of Lemma 16.11 in 
Appendix 6.3, we have 

n<ffo(x,y) ■ $(0,x,y)dxdy 

N 



// ^o(^o,yo)^(^2/)^ /t, ^ 2 ' C/)+ ^ &2 ^ )Ka)d<J ^ 



/ / 0ko( x -9f t >y)fiok(x,y)dxdy+y j / </>{ (x , 27o)$fc(0, ^o, y~o)dx dy . (4.6) 

JO ./<?/< fc=l 

As for the last term of (|4.5p . by solving the backward linear Cauchy problem (14. 2p . we have 

E/' g (T ,i, y) 8 t(T ,o, y> » fr 

^ /-i r9ft ai g f {u f {t-^)) f x 

= > / / 4> k {t - —,l,y)<p ok (x,y)dxdy. 

f^Jo Jo r gf g f 

Since ^ e C$((0, l) 2 ) and t G [0,5] are both arbitrary, we obtain in C°([0, <5]; L 1 ((0, l) 2 )) 
that for k = 1, ■■■ ,iV 



+ 



r aigAufit — £-)) = f T 

— ^ £«_^(t-£.,i, y ), if(:r, y )G[0,^]x[o,l], 

T 9f 9f (4.7) 

. 0fe(*,»>y) = ^ ( X ~9ft,y) = %{t,x,y), if (x,y) e[gft, l]x [0,1]. 



Let us recall the definition of the characteristic £ 3 and also the definition (|4.4j) of (ro,/?o)- 
By solving the backward linear Cauchy problem (|4.2p , noting (|4.7p and (|6.22p of Lemma 16.11 
in Appendix 6.3, we have 



[ I — h-x(-r,l, y)$ k (T,0,y)dydT 
f^Jo Jo 01 



Since <£ G C^((0, l) 2 ) and t G [0,5] are both arbitrary, we obtain in C°([0, 5]; L 1 ((0, l) 2 )) 
that for k = 1 

0, else. 
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For k = 2, • • • , N, we have 



4>t(t,x,y) = < 



2r g/ fc _i(T O ,l ; /3 o ) g _ [X(V 3 .r7)+^(5 3l B / )] if , 
ai9f{u f (T )) 



(4.9) 



0, else. 



For any fixed t G [0, <5], we introduce four new subsets Zj{ , cu 2 ' , u5 3 ' and w 4 ' of [0, l] 2 

-f,t 
u) ± : 

~/,< 
w 2 



5 f t < x < 1, rh(t,5 >t) < y < »fe (*,£/*)}, 
:={(x,i/)| < x <flf/i,0 < y<m(t,x)}u{(x,y)\ g f t < x < 1,0 < y < \ (t, g f t) } , 



w 3'* : ={( x >y)l < X < <?/f, 7?!(t,x) < y < ?? 2 (i,x)j, 
wl' :=[0, l] 2 \(wl' U^' ). 



Here y = r] 1 (t,x) and y = r} 2 (t,x) satisfy 



ds 
ds 



h f(rii,Uf)(s), 7/i(t-— ) = 0, t-^<s<t, 

5/ 5/ 

— 

M 7 ^ «/)(*)> ) = 1 > *- — <S<t- 

5/ 5/ 



For any fixed i G [0, 5] and (x, y) G [0, l] 2 . If (x, y) G w 1 ' , we define £ 5 = (x^,y 5 ) by 

^7=5/, ^ = hf(y 5 ,Uf)(s), £ 5 (t) = (x,y). 

Let us then define (xo,y ) := 1/5 (0))- By solving the backward linear Cauchy problem 

(|4.2p and noting (|6.24|) of Lemma [6. II in Appendix 6.3, the last term of (|4.6p can be rewritten 

as 



iV ,i ,i 



£ 

fc=i 

N 




<fL( x o>yo)®k(0, x , y )dx dy 



o Jo 



fe=i 



I" (x,y) G w 2 ' , we define £ 6 = (x 6 ,y 6 ) by 

~£~ = 5/ 5 ^ = hf(y 6 ,u f )(s), £ 6 (t) = (x,y). 

There exists a unique io such that y^fta) = 0, so that we can define ao := xe(to)- By solving 
the backward linear Cauchy problem (|4.2p . and noting (|6,25p of Lemma [6. II in Appendix 6.3, 
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we have 



N r t r— t 

/ / hf(0,Uf(T))<l> k (T,—X,l)$ k (T,X,0)dxdT 

k=1 Jo Jo «1 

AT „t p _j _ 

■■Y] / / 2 hf(0,uf(t ))(f) k (t , — a , l)$ k (to,ao,0)daodt 

fr[Jo Jo «i 

U) 2 



—f t — _ _ 

If (x,y) G uj 3 ' , we define £ 9 = (x 9 ,y 9 ) by 

~^f = #/' ^ = hf(y 9 ,u f )(s), £ 9 (t) = (x,y). 

There exists a unique tq such that xg(ro) = 0, so that we can define /3 := ^9(^0). By solving 
the backward linear Cauchy problem (j4.2j) and noting (|6.28p of Lemma 16. II in Appendix 6.3, 
we have 



E 

k=2 
N 

E 

fc=2 
N 




JO 



9f(t>k-i(i~, l,y)<f> k (r,0,y)dydT 



9f4>k-i( T o, 1, Po)$k{To, 0, P )d/3 dT 
Jo 



fc=2jV, 



' '3 



Since <£ G C^({0,1) 2 ) and t G [0,5] are both arbitrary, we obtain in C°([0, 5}; L 1 ((0, l) 2 )) 
that for k = 1 



4>i{t,x,y) = < 



r ,_ 

=f/ /T a 2_ 



/o[% 5 >^)+^(y 5l «/)]W^ 



_ 9Sj 



if G Ji\ 



0i(to, — ao,l)e *° 8 " w ' " n ' , if(x,y)Gw 2 

01 

0, else. 



(4.10) 



For k = 2, • • • , iV, we get 



<l>k(t,x,y) = < 



4o^o,Vo)e- J " [ ~ XiV5 ' U)+d -^ ^^ d % if (x,y) G wf*, 

01 (to,— «o,l e Jt <> 1 W8 ' oy w*im) if x,y Gw 2 , 
«i 

lLi(ro, 1 J ) e -4^^ + ^^^)lW^, if (x,y) € 2?, 
0, else. 



(4.11) 



By (|3.1ip . we claim that M{t) G £Is,k and M{t) satisfies the same contraction mapping 
function G. Since M is the unique fixed point of G in £l$Ki therefore we have M = M. 
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Consequently, we have ^ = = 1, • • • , 10), rj i = r]i(i = 1, 2) and rj = 77, so that (xq, = 
(x ,yo), {t ,a ) = (to,ao), (ro,/3 ) = ( r o,A)), = (z = 1,2,3) and lo,{ = ty' (i = 
1,2,3,4). Finally, by comparing the definition (j4T7j) - (p~TT|) of $ f with the definition (I33T1) - 
(|3.25p of (ff, we obtain 0* = (ff. This gives us the uniqueness of the weak solution for small 
time. □ 



5 Proof of the existence of a global solution to the Cauchy 
problem 

Let us now prove the existence of global solution to Cauchy problem (|l.ip - (|1.5p . Noting 
the definition (|3. 1Q[) of Gf and the definition (|3.21j) - (j3.25j) of the local solution <j)f, it is easy 
to check that the following two estimates hold for all t £ [0, 5] 



< M{t) < K, 
,2Ki f aiKx N -f 



2r, 



k V K 2 



T 9f 



(5.1) 

'}• 

(5.2) 

where K is defined by (|3.1|) . K\ is defined by (|3.2|) and K 2 is defined by (|3.3|) . In order to 
obtain a global solution, we suppose that we have solved Cauchy problem (|l.ip - (|1.5p up to 
the moment r £ [0, T] with the weak solution (ftj G C°([0, r]; L p ((0, l) 2 )). It follows from our 
way to construct the weak solution, we know that for any < t < r, the weak solution is 
given by for k = 1 

dh t 



f 



0fo(^o,2/o)e-^ [%2 ' C7)+ ^ fe2 ^ )1(ff)dff , if (x,y) G w{\ 
0, else. 



For k = 2, • • • ,N 



<l>k(t,x,y) :=< 



^ (rr ,2/o)e- / o [%2 ' C/)+ ^ (y2 ^ )1(CT)dff , if G Bp, 

a i9f{ u f{ T o)) 
0, else. 



For = 1 , 



A/ 



b{ (x-g f t,y), if0<i<^-, yG [0,1], 



5/ 



4>{(t,x,y) :-- 



°<Wf(u f (t s_ if £<t<T,y€ [0,1]. 

T 9f 9f 9f 

^o(^o,yo)e- / ° [ ' (y5 ' l/)+ ^ (i/5 ^ )1(<T)dff , if (s,y) G 

<Pi (*o, — «o, l)e 
ai 

0, else. 



(5.3) 



(5.4) 



(5.5) 



■ft [x(ye,u )+1} J-(y 6 ,u f )K*) d *^ z {Xty)€ zf^ (5.6) 
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For k = 2, • • ■ ,N 



Oh , 



4 (x ,y )e-^' u)+ — if (z,y) € wf, 

(/) k (to,—ao,l)e Jt ° l ( ' °y [yb ' fn{ > , (x,y) E u T 2 ' , 

ai _ (5.7) 

^_ 1 (ro,l,/3o)e-^^' C/ ) + ^^^)]W^, if ( x ,y) g £}£*, 

0, else. 



As time increases, for each cell cycle (A; = 1, • • ■ , N) in Phase 1, the characteristic 
passing through the origin may have two possible cases. It may either intersect with the 
front face (see Fig 6 (a)) or intersect with the right face (see Fig 6 (b)). For each cell cycle 
k (k = 1, • • • ,N) in Phase 3, the characteristic passing through the origin will definitely 
intersect with the front face due to the fact that hf(l,Uf) < (see Fig 6 (c)). We get in any 
case that 

N f f ~f N f 1 f 1 

J J k=1 Jo J0 



k=l 
N 



2 



N N 

+}2 / / a 2 jo(joy-i^s)(f>{(t,x,y)dxdy+^2 / / a 2 jo('joy-+ls)<t>{(t, x,y)dxdy 
k=\t/ >t k=\tl t 



( t ( 




(a) (b) (c) 

Figure 6: For t large enough. Case (a): the characteristic passing through the origin intersects 
the front face at (to, 1,/3o) in Phase 1; Case (b): the characteristic passing through the origin 
intersects the right face at (to,ao,l) in Phase 1; Case (c): due to the fact that hf(l,Uf) < 
in Phase 3, the characteristic passing through the origin will definitely intersect the front face 
at (i,l,/3 ). 

t 

We can prove that estimate (|5.ip holds for every t € [0, T] by tracing back 4> k (t,x,y), 
<p k (t,x,y) and (f> k (t,x,y) along the characteristics to the initial data at most N times. More- 
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over, noting the definition (|5.3|) - ([5.7p of the giobai sohition, it is easy to check that the uniform 
a priori estimate (|5.2p holds for every t £ [0, T\. Hence we can choose 5 £ [0, T] independent of 
r. Applying Lemma 13.21 and Lemma |4 . 1 1 again . the weak solution <f>f € C°([0, r]; L p ((0, l) 2 )) 
is extended to the time interval [r, r + 5] f] [t, T]. Step by step, we finally obtain a unique 
global weak solution <f>t € C ([0, T]; L p ((0, l) 2 )). This finishes the proof of the existence of a 
global solution to Cauchy problem (|1.1|) - (|1.5|) . 



6 Appendix 

6.1 Introduction of the model 

In this section, we first recall the systems of equations describing the dynamics of the cell 
density of a follicle in the model of F. Clement [H [6]. The cell population in a follicle / is 
represented by cell density functions k (t,a,j) defined on each cellular phase Qj k with age 
a and maturity 7, which satisfy the following conservation laws 

dcj>[ k d(g f (u f )<l> f jk ) a(h/(7,"/)4>j, fc ) rn ,/ . _ 

-QT + da + 57 = - X b>m# ™QU (6-1) 



Qj,k := x [°> T]-> for J = 1, 2, 3, / = 1, • • • , n, where 

:= [(k - l)o 2 , (A; - l)a 2 + ai] x [0, 7J, fi 2 ,fc := [(*; - l)o 2 + 01, A;a 2 ] x [0, 7J, 
^3,fe := [(A; - l)a 2 ,fea 2 ]. 

Here k = 1, • • ■ , iV, and iV is the number of consecutive cell cycles (see Fig 7). 
Define the maturity operator M as 

M(tp)(t):= / ~ /i p(t,a n )dad 1 . (6.2) 

Jo Jo 

Then 

3 N 



M f :=J2J2 / ld k (t,a, 1 )dad 1 (6.3) 

T^t~,Jo Jo 



j=i k=i ' 

is the global follicular maturity on the follicular scale, while 



3 N 

' Jmax I u )7m.r 



M:= EEE/ / ^ jk (t,an)dach (6.4) 

is the global maturity on the ovarian scale. 

The velocity and source terms are all smooth functions of uj and U, where Uf is the local 
control and U is the global control. In this paper, we consider close/open loop problem, that 
is to say Uf and U are functions of Mf, M and t. As an instance of close loop problem, the 
velocity and source terms are given by the following expressions in [1] (all parameters are 
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positive constants) 

gf(uf)=T g f(l-g 1 (l-Uf)), in Q 1)k , 
g f (u f ) = T gf , in % jfc , j = 2,3, 

— U j 

hf(qf,Uf) = T hf (--f 2 + (ci7 + c 2 )(l - e~*~)), in % ifc , j = l,3, 

hf(j,Uf) = A(7,J7) = 0, in S7 2ifc , 

A( 7 ,C/) = Ke- (2 ^ )2 (l-f/), in i = 1, 3. 

J>2Mf I — U 

u f = b(M f )U = mm{ bl + -£-,!}■ Po + U s + 1 + eC(M l m) ]- 
The initial conditions are given as follows 

^ fe (0,a,7) = ^o(^7)k fc , J = 1,2,3. (6.5) 
The boundary conditions are given as follows 

f J 2T g f4>2 fc-i(*' ~~ l) a 2, 7), for fe > 2, 

5 f /(«/)</ , i/ c (^(^-l)a2,7) = { ' for 7 G [0,7 fl ]. 

[0, for k = 1, 

0{ fc (t, a, 0) = 0, for a G [(fc - l)a 2 , (A; - l)a 2 + ai], 

T s/^2,fc(*' ( fc ~ 1 )°2 + a i,7) = 9f(u f )(f){ k (t, (k - l)a 2 + 01,7), for 76 [0,7 S ], 

02 a, 0) = 0, for a G [(k — l)a 2 + ai, fca 2 ], 

r f — 1)02,7), for A; > 2, 

03fc(*' ( k ~ 1 ) a 2,l) = { ' for 7 e [7s, 7m], 

I 0, for k = 1, 



0{ fc (t, a, 7 S ) 



0{ )fc (t,a,7 s ), for a G [(/c - l)a 2 , (fc - l)a 2 + 01], 
0, for a G [(k — l)a 2 + 01, fea 2 ]- 



For sake of simplicity, we denote by 70 := 7 m — 7s > in the whole paper the difference 
between the maximum 7 S of the maturity in Phase 1 and the maximum j m of the maturity 
in Phase 3. 

6.2 Mathematical reformulation 

In this section, we perform a mathematical reformulation of the original model introduced 
in Appendix 6.1. Notice that in system (|6.ip . all the unknowns are defined on different 
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' . 4^2 


< 

i . . . 


' . 4, k 


< 

* . . . 


' „ 4,N 


4,1 


' 4,i 


- 4,2 


< 




4, k 


' 4,k 




4,N 


t 

4,N 



Q o.\ a 2 a 2 + fli 2a 2 (fc — l)a 2 (fc — l)a 2 + a^a 2 (iV — l)a^Af — l)a 2 + ajVa 2 



Figure 7: The illustration of iV cell cycles for follicle /; a represents the age of the cell and 7 
represents the maturity of the cell. The top of the domain corresponds to the differentiation 
phase and the bottom to the the proliferation phase. 



domains, so that one has to solve the equations successively. Here we transform system (|6. 1 j) 
into a regular one where the unknowns are defined on the same domain [0,T] x [0, l] 2 . We 
denote by 4 the density functions and by g^ and hf the age and maturity velocities for Phase 
1; 4>l the density functions and gf the age velocities for Phase 2; <f>l the density functions 
and gf,hj the age and maturity velocities for Phase 3. 
Let 

</>*:(*> y) := <j>{ )k {t, a, 7), (a, 7) G 9, 1>k , (6.6) 

where 

X := a - {k - 1)a \ y:=^ (6.7) 
ai 7 S 

so that we get (0,7) G Q,^ k (x,y) G [0, l] 2 , and 

($fe)t + (9f4k)x + (hf4k)y = -A (6.8) 

with 

?/(«/):= V(y,«/):= M7 ' y,U/) , A( y ,C/):=A( 7sy ,C/). (6.9) 

J ai 7 S 

Let 

y) : = 02,fe(* 5 a , 7), («, 7) G ^2,fc, (6.10) 

where 

a - (fc - l)a 2 - ai 7 
x := , y := — . (6.11) 

a-2 ~ ai 7 S 

So that we get (0,7) G n 2 & G [0, l] 2 , and 



a/v -i-?;.(2f\ — n ^. — IkZ 

Let 



(<^)t + gf{4) x = 0, (6.12) 

«2 — a l 



4 k (t,x,y) ■= ^3, fc (M>7), («, 7) G ^3,fc, (6.13) 

where 

x:= a - (fc - 1)a2 ,,:=^, (6.14) 
02 7o 

so that we get (0,7) G f^fc -4=> (a;, y) G [0, l] 2 , and 

+ (9f4 k h + (h f 4 k )y = -W k , (6-15) 
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with 

9f •■= —, hf(y,Uf) := — — , \{y, U) := \{j y + -y s , U). (6.16) 

0,2 70 



Accordingly, let us denote the initial conditions (|6. 5[) with new notations by 

$fo = (4>io( x iy)r-- ,~<i>No( x iy), <i>io( x ,y)r ■ ^L( x ^y)r-- ,4>N ( x ,y)), (6.17) 

where 

fao( x >y) : = W,fc(°; a ,7), <i>ko( x >y) := ^2,fc(°' a '7), := 4,fe(°' a '7)- (6-18) 
Let 0/ = (^!,-- - ,^at ; <£{>••• ,<£at> ^fr- - ><^at) T > w e have 



<M*> ^ y) t + (Af<f>f(t, x, y)) x + (Bf<j> f (t, x, y)) y = C(p f (t, x, y), (6.19) 
te[0,T], (x, y) G [0, l] 2 , 



where 



Ar N N 



A f : = diag {g f , ■ ■ ■ ,g f , g f ,-- - ,g f g f , ■ ■ ■ ,g f }, 

N N N 

B f : = diag {hf,--- ,hf, 0, • • • ,0 hj,--- ,hf}, 



N N N 



C : = -diag {A, • • • , A, 0, • • • ,0 A, • • • , A}. 
From the original expression of hf(^,Uf) (see Appendix 6.1), let us define 



. , ci(l-e s )±Vq(l-e « ) 2 +4c 2 (l-e s ) 
7±(«/) : = j ' 

It is easy to see that 7+ (it/) is an increasing function of Uf, and 

hf(rf,u f ) = T hf {~i + {uf) - 7)(7~7-(«/))- 
Hence, when 7 = 7 + (n/), we have /i/(7, «/) = 0, moreover 

hf(j,Uf)>0, if < 7 < 7+C"/), 
h f (j,Uf)<0, if7>7+(u/). 

Furthermore, under the assumption that the local control iiy satisfies 7+(itf) > 7 S , and 



< 7+0/) < 7+(l) ) < 7m- 

From the fact that < 7 < 7 S in Phase 1, and 7 S < 7 < 7 m in Phase 3, we have 

7i/(7,u/)>0, 0<7<7 S , 
hf(%,Uf)>0, hf(j m ,Uf) < 0. 
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Corresponding to the new notations, we have assumptions 



h f (y, Uf )>Q, VyG[0,l], 

hf(0,u f )>0, hf(l,u f )<0. (6.20) 



Since 



M (W,fc) = // 74>lk(t,a,7)dad>y = J au 2 y(p{(t,x,y) dx dy, 



M (^2fc) =// 702 fc^ 13 , l)dad-f = I [ (a 2 - a 1 )j 2 y4>{(t,x,y)dxdy, 



o j o 



2.1: 



M ( ( l ) {k) = J J 7${,k( t ' a >i) dad 'y = J J o a 27o(7oy + J s )<f> f k (t,x, y)dxdy. 

We have 

N N N 

M f (t) = M(< fc ) + ]T M(4 ik ) + M(4,k) 
k=l k=l k=l 

N ,1 ,1 . JV ,1 ,1 



/ / ai^y4>{(t,x,y)dxdy + J2 [ [ (o 2 - a 1 )j 2 y4> f k (t,x,y)dxdy 
k=i Jo Jo k=i Jo Jo 

+ 52/ / a 2lo(loy + ls)4>{(t,x,y)dxdy. 
r - , jo jo 



fc=1 J0 jo 

6.3 Basic lemmas and some notations 

The following lemma is used to prove the existence of the weak solution to Cauchy problem 
(|l.ljl - (jl.5p . when we derive the contraction mapping function G, change variables in certain 
integrals (see Section 3) and prove the uniqueness of the solution (see Section 4). 

Lemma 6.1. The characteristic £2 = (2:2,2/2) passing through (t,x,y) intersects the bottom 
face at (0, xo,2/o)- We have 

d( X >y) _ Jo ^f{y2,u f )(a) da /g 2] x 

d(x ,y Q ) 

The characteristic £3 = (2:3,2/3) passing through (t,x,y) intersects the back face at (tq,0,Pq). 
We have 

= -g f (u f (r )) ■ A ^<"'«/>W* (g.22) 

The characteristic £4 = (2:4,2/4) passing through (tq, 1,/?o) intersects the bottom face at 
(0, xo, J/o) • We /iowe 

5(a:o,yo) g f (uf(r )) 
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The characteristic £5 = (3:5, 2/5) passing through (t, x, y) intersects the bottom face at (0, xq, i/q). 
We have 

9 ( X >y) = Jo ^f(ys,u f ) (<t) da , g 
5(^0,2/0) 

T/ie characteristic = (xq,i/q) passing through (t,x,y) intersects the left face at (to,otQ,0). 
We have 

o{to,ao) 

The characteristic £7 = (2:7,2/7) passing through (to,cto,l) intersects the bottom face at 
(0, XQ,yo). We have 

d(x Q ,y ) hf(u f (t ),l) 

The characteristic = 0^8 1 2/8 ) passing though (to,ao, 1) intersects the back face at (tq,0,/3q). 
We have 

d(ao,t ) _ -g f (u f (ro)) j<o a ^± {ys , Uf){a)da 

d{r ,p ) h f (u f (t ),l) 
The characteristic £9 = (2:9,2/9) passing through (t,x,y) intersects the back face at (ro,0,/3o). 
We have 



g fe lrt>SirW a 'W r > da . (6.28) 

o(t ,(3o) 

The characteristic £10 = (2:10,2/10) passing through (tq,1,(]q) intersects the bottom face at 
(0, XQ,yo). We have 

d(x ,yo) gf 

The proof of this lemma is trivial, and can be found in [TJ] of 1-space dimension case, we 
omit here. 

The expressions of contraction mapping coefficients C( and c£ are as following 
c{ : ^ m{ 2Kl ~ «^±«3* + ^ £ (0.30) 



k=l 



2(a 2 -a l )-,l(iq + 2tKlK 2 +iK l K 2 )^ l ~, a 27 g,(2K 1 + 2tg 1 2 )^ nIJ 



ZV i; fc=2 1 k=l 



/ _ ai 7 4(2^ + 2Ki + 12^x^2 - 2ti^if 2 ) ^ „-/ 

G 2 - jr,d-*iro 



k=l 

2(a 2 - ai)7^.(^ + §K X K 2 ) ^ ~ 4a 2 7^^i M r/ 



+ " ^(l-t^O ^2 1-^1 
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